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The total Born approximation cross section for pair production in a coulomb field given originally 
by Racah is used to derive two simple and rapidly convergent analytic expansions for this cross sec- 
tion, one valid for high energies, A:, of the incident photon, the other for energies near the threshold. 
For k >4<mc 2 the fractional error involved in using the first four terms in the high energy expansion is 
<4.4X 10" 5 . For k < 4 rnc 2 the fractional error committed in using the first five terms of the low energy 
expansion is < 1.1 X 10 4 . The leading terms in these expansions are the well known high energy and 
low energy limits of this cross section, respectively. 
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The cross section for the production of an electron-positron pair in the field of an unscreened 
Coulomb potential V—~Ze l \r was first given in first Born approximation in a number of papers 
nearly thirty five years ago by W. Heitler and F. Sauter [I], 1 H. Bethe and W. Heitler [2, 3], and 
G. Racah [4J. References 1, 2, and 3 give this cross section, differential in energy but integrated 
over all particle directions, without any approximation in the energy [5] as well as in a form valid 
when the energy of each of the pair particles is extremely relativistic (E.R.) [6], (d > 1, €z > 1): 
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^here 0=^- (— J = Z 2 ■ 5.7938 X 10~ 2H cm 2 . (2) 

and 6i, €2, k are the energies of the positron, electron and photon respectively, in units of mc 1 . 
References 2 and 3 give the total cross section (eq (1) integrated over d), but only in the extreme 
relativistic limit (k > 1): 
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Only in the papers by Racah can one find an expression for the total cross section involving no 
high energy approximations. Unfortunately, this expression involves indefinite integrals over 
elliptic integrals, so that it is quite inconvenient for obtaining numerical values of the total cross 
section. In fact, for this purpose it is equally simple to integrate directly, numerically, the differ- 
ential cross section [5] in which no high energy approximations have been made. However, both 



1 Figures in brackets indicate the literature references at the end of this paper. 
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for numerical evaluation and for an understanding of the analytical approximations involved in 
the extreme relativistic expression (3), it is clear that it would be convenient to have a simple 
analytic expression for the total cross section, valid over a large range of energies, which reduces 
in an obvious manner to (3) in the limit of large A. A step in this direction was made in a paper by 
Hough [7] who gives a corrected expression for the total cross section: 
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to be applied at lower energies than those for which (3) may be used. In addition, a number of 
empirical analytical expressions are given to cover the region 2 ^ k ^ 15. They are, however, not 
really very simple, and they do not go smoothly into the high energy expression. Indeed, while 
the expression (4) is, as far as numerical values are concerned, an improvement over (3) for mod- 
erate energies, it is rather deceptive from the analytical standpoint: The corrections to (3) for 
moderate k are not of the form \jk. The expression (4) results in fact from an integration of the 
high energy differential cross section (1) between the limits €=1 and e 1 = k— 1 rather than 
and A, which leads to (3): 
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in which 12-8 In 2 = 6.45 . . . 

However, the extreme relativistic differential cross section in the integrands here has errors of 

relative order 1/cf and 1/cf, and these may be expected to give contributions of order 



l + f € l€2 



A- 3 1 



3 ' 7 /, 2€i€ 2 1\ , „ /lnA\ „ (I 



i.e., of exactly the same order as the correction calculated by changing the limits of integration 
from (0, k) to (1, k— 1). This is in fact borne out by the expansion given here, in which we see that 
no terms of the form 1/A actually appear in successive terms in the high energy expansion. 

In this note we return therefore to the expression for the total cross section without any ap- 
proximation in the energy, given by Racah, and derive from it two relatively simple and rapidly 
convergent expansions — one for large A, the other for k near threshold (k = 2). The first few terms 
in each of these expansions are shown to suffice for both analytic and numerical purposes for all A. 

The results which we obtain are the following: For A near threshold we derive in fact two 
series, one in the expansion parameter 
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the other in the somewhat more complicated expansion parameter 
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This second series has the advantage of being more rapidly convergent. These expansions are 
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and 
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For large k we derive an expansion in the parameter 
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The convergence of this high energy expansion is extremely rapid, even for rather low energies. 

For example, for k = 4>, the fractional error, 6, (cr a p Pr()X = cr exact (l + 8)) involved in using the terms 

shown explicitly in (12) is only + 4.4 X 10~ 5 . The fractional errors obtained using the terms shown 

explicitly in (9) and (10), also for k = 4, are — 1.1 X 10 :i and — 1.1 X 10~ 4 , respectively. For k > 4, 

the errors in (12) become smaller and for k < 4 the errors in (9) and (10) become smaller. 

The first term in each of these expansions has been given previously by Racah [8]. Their 

regions of validity are shown in figure 1, where we give graphs of the ratio of the first term in (10) 
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total cross section, cr. It is to be noted from figure 1 that the first term of (10), valid for k near 
threshold, is more than 10 percent ofT for k > 1.5 MeV, and that the first term of (12), valid for 
k very large, is more than 10 percent off for k < 10 MeV. In figure 2 we give a graph of the total cross 
section, cr, in millibarns, as a function of the incident photon energy, /r, in MeV. In table 1, values, 
(cr/(/)), of the total cross section in units of (/) are given for photon energies ranging from 1.03 to 
1000 MeV. 

We proceed now to the derivation of these expansions, our starting point being the exact 
Born approximation cross section as given by Racah: 
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FIGURE 1. The ratio of the first term in expansion (10) ( C r t * reshold = </> • — ( j ) 

/ — /28 218W 

and of the first term in expansion {12) fcrf-*- = <f) I — In 2k — —J J, to the exact tota< 

cross section, cr, as a function of the incident photon energy k, in MeV. 
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FIGURE 2. The total pair production cross section, cr, in millibarns, as a function of 
the incident photon energy, k, in MeV. 



where K{j]) and E(r\) are the complete elliptic integrals of the first and second kind, respectively 
[9J, and are denned by 
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TABLE 1. The total cross section in units of cf> as a function of incident 
photon energy in MeV 



k (MeV) 


<r/<J) 


k (MeV) 


o-/c/) 


1.03 


9.788 X10- 7 


10 


3.647 


1.05 


3.994 X 10- 5 


15 


4.782 


1.10 


7.607 X 10" 4 


20 


5.618 


1.15 


2.978 X 10 ; < 


30 


6.826 


1.20 


7.134 XlO" 3 


40 


7.696 


1.25 


1.342 X10- 2 


50 


8.377 


1.30 


2.189X10" 2 


60 


8.936 


1.40 


4.507 X 10- 2 


80 


9.822 


1.50 


7.577 X10-* 


100 


10.511 


1.75 


1.777 X10" 1 


150 


11.767 


2.00 


3.030X10-' 


200 


12.659 


2.50 


5.836x10' 


300 


13.919 


3.00 


8.716X10-' 


400 


14.813 


4.00 


1.416 


500 


15.507 


5.00 


1.9O0 


600 


16.074 


6.00 


2.328 


800 


16.968 


8.00 


3.052 


1000 


17.662 
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where 
and 



(a) n = a(a + 1) . . . (a-\-n-\) for n > 
(a) =L 



Thus for r] < 1, we may write out the power series for K and E in (13), integrate term by term, and 
collect successive powers of 7). This straightforward procedure leads to our expansion (9). 
Expansion (10) is obtained by using the transformations [9] 
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The series for K(k) and E(k) converge faster than those for K(r)) and E{iq) since X < tj. Sub- 
stituting (16), (17) and (18) in (13), we write out the power series for K(k) , E(k) and tj (a series in 
X 1 ' 2 ), perform the integrations term by term with respect to X, collect successive powers of X 1 2 
and write the series in terms of the variable p = 2k li2 = 2rjl[l + Vl — -n 2 ] since p^r) as tj — » 0. 
This gives the expansion (10). 
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In order to obtain the high energy expansion (12), we note that 17—1 for k> 1. We therefore 
need expansions of the functions K{rj) and E(r)) about 77 = 1. To this end we rewrite the trans- 
formations (16) and (17) with X and r\ interchanged: 
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and write (19) and (20) in the form 
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In (13), the terms preceding the integrals can now be rewritten conveniently in the form 
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The desired expansions of K and E are 
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where i//(z) = r'(z)/r(z), and, in particular, [10) 
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From (28) we may write 
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Equations (26) and (27) are then written in the form 
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Substitution of (30) and (31) in the right hand side of (25) gives the desired high energy expansion 
of the terms in (13) preceding the integrals. We note that the terms in this expansion are either 
of the form of a power series in /x 2 = (2/A) 2 or In 2k times a power series in fx 1 . 

We next consider the high energy expansion of the integrals in (13). Each of these integrals 
has the factor 
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Now from (19) we have 
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Since we want an expansion of /t for k> 1, i.e., 2/k< 1, we separate that part of the integrand 
which becomes singular as z — > 0. From (30) 
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as z^O. Thus we may write 
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The evaluation of the definite integral in (38) is given in the appendix. The result is 



f 1 dz 

Jo z 



£(Vl-z 2 )-ln 



2 In 2 2-^tt 2 . 



(39) 



In the last integral in (38) we substitute (30), with fx replaced by z, 2k by 4/z, and integrate term by 
term. We then obtain 
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Finally, we consider the second integral in (13): 
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upon integrating by parts. Again we substitute (33) and (35). and obtain 
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As with /r, we write 
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The evaluation of the first integral here is also given in the appendix. The result is 
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In the last integral in (43) in we again substitute (30) and integrate term by term. Since these in- 
tegrals appear in (13) only in the form of the sum /] + / 2 , we give this result: 
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Appendix 

We evaluate here the definite integrals appearing in (38) and (43): 
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with a 3= 0, so that 
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Substituting £ = 1 — z 2 in (A.3) and using (14b), we have 
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Here 
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We then expand this expression in a power series in a, keeping, in the final result, terms linear 
in a. We note 

r'(z)=««z)r( Z ) 

r"(z) = xl,'(z)nz) + iljHz)r(z) (A.9) 

and the following values for specific values of z: [10] 
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We find 
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from which J\ and Jz follow from (A.4) and (A. 5). 
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We would like to express our thanks to J. H. Hubbell for performing the computer calculations 
involved in the preparation of figures 1 and 2 and table 1. 
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